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Abstract 

■ In this paper we investigate the power spectrum of thermal fluctuations in very early stage of 

(N 

5h ■ viscous cosmology. When the state parameter as well as the viscous coefficient of a barotropic 

o3 

fluid is properly chosen, a scale invariant spectrum with large non-Gaussianity can be obtained. 
In contrast to the results previously obtained in string gas cosmology and holographic cosmology, 
we find the non-Gaussianity in this context can be fc-independent such that it is not suppressed at 

'o'. 

O" 1 large scale, which is expected to be testified in future observation. 
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I. INTRODUCTION 



Recent observed Cosmic Microwave Background (CMB) anisotropy[l-3| and the large 



scale structure of the universe 



4j-|8] can be viewed as the results of the primordial den- 



sity perturbations, which is characterized by a nearly scale invariant and Gaussian power 
spectrum [9]. Though at present completely understanding the origin of the primordial den- 
sity perturbations is still an open question, one traditional point of view is thinking of 
quantum fluctuations in vacuum as the seed of such classical perturbations. Nevertheless, 
there is an alternative conjecture arguing that it may be due to thermal fluctuations of 
matter sources during the inflationary stage 10Ml3l|. Unfortunately, thermal fluctuations in 
standard inflation models always generate a power spectrum with spectral index n s = 4 
{14! [fo] |. Recently progress has been made to overcome this difficulty People show that 
once some new physics is introduced, for instance in the context of the noncommutative 
inflation, holographic cosmology and loop cosmolo gy, the scale invariant spectrum can be 



implemented in the thermal scenario as well 
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20]. 



Recently the thermal origin of the primordial density perturbation has received more 
attention since the non-Gaussianity of CMB has been further disclosed in the latest ob- 
servation data. It indicates that at 95% confidence level, the primordial non-Gaussianity 
parameters for the local and equilateral models are in the region —9 < f^l al < 111 and 



— 151 < fx™ 1 < 253, respectively [9]. If this result is confirmed by future experiments 
such as the Planck satellite, then it will be a great challenge to many slow-roll inflation 
models since the Non-Gaussianity has to be greatly suppressed at |/jvz| < 1 in most of 



these models 



2ll . [22| . Contrasting to fluctuations originated from vacuum, the fluctuations 
in thermal scenario are always not strictly Gaussian. This may provide an alternative way 
to seek the observable non-Gaussianity. Therefore, the thermal non-Gaussianity has been 



widely discussed in recent literatures 
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In this paper we intend to investigate the power spectrum of thermal fluctuations in 
the very early stage of viscous cosmology, where the matter source is a viscous barotropic 
fluid (the state parameter w = - £ [— 1, 1] is a constant). Since the viscous effect has an 
impact on the perturbation modes after they cross the thermal horizon, it is expected that 
the power spectrum should be different from the previous results even without other new 
mechanism or extra structure introduced. As a matter of fact, we find that provided the 



2 



state parameter as well as the viscous coefficient of a barotropic fluid is properly chosen, 
a scale invariant spectrum with large Non-Gaussianity can be obtained in this framework 
indeed. 

The outline of our paper is the following. In next section we present a brief review on 
the thermal fluctuations in the FRW universe filled with a perfect fluid and demonstrate the 
difficulty of gaining a scale-invariant spectrum in this model. Then in section III we show 
that this difficulty can be overcome by introducing a viscous fluid and properly choosing its 
viscous coefficient. Its non-Gaussianity is investigated in section IV. Finally the comparison 
with other cosmological models with thermal fluctuations is given in the section of Conclusion 
and Discussion. 

II. A REVIEW ON THERMAL FLUCTUATIONS IN THE FRW UNIVERSE 
WITH A PERFECT FLUID 

In this section, we briefly review the thermal fluctuations in very early universe and derive 
its power spectrum based on thermodynamical consideration. First of all, we consider the 
FRW universe filled with a perfect fluid with state parameter w. For a perfect fluid, the 
thermal system can reach equilibrium state through interactions at sub-Horizon scales. In 
the thermal scenario of cosmological fluctuations, this process is called thermalization. Since 
during this process both of the extensive internal energy U and entropy S are well defined 
for a global equilibrium state, we have the following form of the first law of thermodynamics 



For barotropic fluid, one finds that the relation between the energy density and the temper- 
ature is uniquely fixed by this thermodynamical lawfl^. [loj 



where m = 1 + — and A is an integral constant. For radiation with w = 1/3, this equation 
is nothing but the famous Stefan-Boltzmann law. On the other hand, the partition function 
of the thermal fluid is defined as 



dU = TdS - pdV. 



(1) 




p = AT m , 



(2) 




(3) 



r 
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where /3 — T 1 . The internal energy U inside a volume V is given by 

Therefore, the two-point correlation function for the density fluctuation Sp = p — (p) can be 



obtained as 



(jg) _ (£ 2 ) - (E) 2 1 tZ 2 logZ J_dC/ 
V 2 y 2 ~ 



T 2 C V 
R 6 



(5) 



where CV = (f^)v = V'^f = V^p' is the specific heat and R ~ W is the size of the thermal 
horizon. Thermal fluctuations generated from the matter inside R can be described by the 
thermodynamics above. But when thermal modes are pushed outside the horizon, they are 
frozen and become non-thermal governed by the theory of perturbations. Next we calculate 

iteratures we identify the thermal 



the power spectrum of perturbations. Followin g m ost of 
horizon R with the Hubble horizon H 
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271 ]. which means our calculation 



of spectrum is always taken at the Hubble scale. 

If perturbations are deeply in the horizon, the 0-0 component of the perturbative Einstein 
equation will reduce to the Poisson equation which relates the curvature fluctuations and 
the density perturbations Spk as [35] 



/c 2 $ fc = 47rGa 2 5pk, 



(6) 



3 

where 5pk = k~?5p. Thus the power spectrum can be obtained by combining Eq.([3]) and 



k 3 



(**) ~ T T2 P' 



(7) 



2tt 

where we have used the condition R = H~ l = |. When the $ modes leave the horizon, 
their amplitudes get fixed at whatever thermal amplitudes they have at crossing k = aH. 
For simplicity, we consider the spatially flat universe. Then using the Friedmann equation 
if 2 oc p, we have 



V*(k) 



T 2 p> 
s/P 



(8) 



k=aH 



For a constant w, substituting the thermodynamical relation fl2]) into above equation leads 
to 

dlnV, 1 + m (g) 



cilnT 



m 
~2' 
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Furthermore, from the conservation equation of the fluid one has p oc a ^ l+w ) such that 



a oc T 3 ( 1+m > 



(10) 



Thus we have 



d In k m(l + 3w) 



(11) 



d\nT 6(1 + w) ' 

where we have used k = aH oc a^fp. Using Qj and fTTTj) . it is straightforward to calculate 
the spectral index as 



ng-1 



d In 'P.j, d In 'P.j, <iln T 2 + m w; + 1 

3 ■ 



(12) 



dliak dlnT dink m 3w; + 1 

In this equation the relation m = 1 + — has been applied. Therefore th e sp ectrum is always 



blue and independent of the value of w. Of course this ll no-go result" 15| is not consistent 
with the current experiments, in which n s is restricted at n s = 0.9601qo}3 js]. Then, to 
obtain a scale invariant spectrum in thermal scenario, one need introduce new physics to 
either relax some constraints due to thermodynamics, for instance as presented in 14|, Il9| . 
or modify the standard cosmological equations, as presented in Different from above 

considerations, in next section we would like to argue that the no-go result above can also 
be avoided if the matter source is a viscous fluid rather than a perfect one. 



III. REALIZATION OF SCALE-INVARIANT FLUCTUATIONS IN VISCOUS 
COSMOLOGY 



Viscous cosm ology h as been widely applied to investi gate the structure and the evolution 



of the universe 
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-|50| (For recent review we refer to 5l|)- In this context the energy 



momentum tensor of the fluid is given by 



T„ v = (p + V ~ ^H)u^u u + (p - 3£H)g^ lu , 



(13) 



where the bulk viscosity coefficient £ = £(p) is usually a function of the energy density 
p of the fluid. Moreover, as pointed out in 52[, £ should be positive if the second law 
of thermodynamics is respected. In FRW universe the conservation equation V^T^j, = 
becomes 



p + 3H(p + p-3£H) =0. 



(14) 
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In this paper, we choose a special kind of viscous fluids with a bul 



t viscosity coefficient 
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53j. In this case, 



oc pa, which has already been investigated in some literatures 

= a is a positive constant which is relevant to the viscosity of the fluid. Then given a 
constant w, we have a relation between the energy density p and the scale factor from the 
equation above 

p oc a- 3{1+w - a \ (15) 

Obviously, the evolution of the universe depends not only on the state parameter w but also 
on a. 

Now we discuss the thermalization of the matter source. Usually, the thermal scenario 
of cosmological fluctuations is established on the thermodynamics of the equilibrium state. 
Unlike the case with a perfect fluid, the viscosity brings dissipative effect such that the 
global description of the first law of thermodynamics in Eq.flT]) no longer stands in general. 
Nevertheless, each particle in a viscous fluid satisfies the Gibbs relation in a local equilibrium 
state 

Tds = d^+pd-, (16) 
n n 

where s is the entropy and n is the particle number in a local equilibrium element. For the 
barotropic case, it can be rewritten as 



ds = ^dp - 9 * + W dn. 



(17) 



Employing the integrability condition 

-(--) 

T \dnn/ p 



w 



n- 



d p 

dpfJn 



(18) 



one can still derive the Stefan-Boltzmann law p oc T m with m = 1 + 1 [45|, |54|. This result 
indicates the viscosity does not change the relation between the energy density and the 
temperature of the fluid. 

In a heuristic manner we may argue that the anisotropy and the inhomogeneity of our 
universe totally originate from the fluctuations in such a local equilibrium element. Consider 
a quanta of the fluctuation with a momentum P r . Since the momentum uncertainty of a 
quantum particle AP r is of order of P r 55m7l|. from the uncertainty principle we have 

h 



Pf — f%Jv> 



rph 



AP r ~ 



Ax r 



(19) 
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where k rp h — ^ is the physical wavenumber. Then we have Ax r ~ j*-. Using the crossing 
condition k r = a r H r , we derive its position uncertainty Ax r ~ H~ x . Since all of the 
fluctuations happen in the local equilibrium element, the position uncertainty should not be 
larger than the scale of the element Aa; r < R. We find 

H~ l ~ Ax < R, (20) 

when each k crosses the horizon. In other words, this local equilibrium element within a 
radius of the Hubble scale can be viewed as the cradle of the cosmological perturbations. 
Moreover, we point out that both the background equation if 2 oc p and the Poisson equation 
OH]) will not change due to the viscosity effects and this can be understood from the fact that 
the viscous term does not appear in the 0-0 component of the energy-momentum tensor, 
as shown in Eq. (ll3p . Therefore the correlation functions can be similarly calculated as in 
previous section. Inserting the relation in Eq.(j2]) into ( TT5|) . we find the variation of the scale 
factor with the temperature is modified as 

— m 

oocT3(i+«). (21) 

Repeating the calculation in previous section, we derive 

dhik m{l + 3w — 3a) 
d\nT 6(1 + w — a) 

As a result, the spectral index for a viscous fluid can be obtained as 



(22) 



dlnV$ 2 + m w — a + 1 3w + 1 w — a + 1 

n s — 1 = — — — = 3 = 3 . (23) 

dink m 3w — 3a + 1 w + 1 3w — 3a + 1 

We find that the effect of the viscosity does imprint on the power spectrum of the cosmolog- 
ical perturbations. Especially, if we choose w ~ a — 1 or w ~ - 1 with any non- vanishing a, 
a nearly scale-invariant power spectrum of the perturbations can be implemented. The case 
with w ~ a — 1 is nothing but the standard inflation, because H = — 47rGp(l + w — a) ~ 
implies a nearly de Sitter expansion. When w « — |, since a > as required by the second 
law of thermodynamics, we have - = — 4irGp(l + 3w — 3a) > 0. In this case, the universe 
undergoes an accelerating expansion as well. Therefore, in contrast to models with a perfect 
fluid, the thermal scenario in viscous cosmology may provide an alternative interpretation 
on the scale invariance of the power spectrum. However, to be well consistent with the 
current observation constraints, we need further study the non-Gaussianity in this thermal 
scenario, and this is what we intend to do in next section. 



IV. THERMAL NON-GAUSSIANITY IN VISCOUS COSMOLOGY 



Another important quantity which is expected to be measured precisely in coming ex- 
periments is the non-Gaussianity of CMB. The amount of the non-Gaussianity is usually 
estimated by the quantity /jvl which is defined by the three-point correlation function for 
the curvature perturbations (Ck 1 Ck 2 Ck 3 )- Generally, two cases are considered, one is the 
local non-Gaussianity labeled by f^f? 1 with k\ <C &2 ~ ^3? an d the other is the equilateral 
non-Gaussianity denoted by ft?T l with k\ ~ k<i ps k^. More detailed discussions on these 



quantities can be found in 



58 



59[. Recent WMAP-5 data shows the possibility that there 



exist large non-Gaussianities of cosmological perturbations [9[. If this is true, thermal sce- 
nario may give a more consistent result. In this section, we will calculate the equilateral 
non-Gaussianity of thermal fluctuations in the context of viscous cosmology. 

The non-Gaussianity is quantificationally described by the non-Gaussianity estimator 

requil _ 5 3 ((f) 

tNL ~ is" (CI) (CI) ' (24) 

where Cfc is the curvature perturbation on slices of uniform energy density, while (Cf ) and 
(C|) are two-point and three-point correlation functions, respectively. Since no entropic per- 



601 ]. the cosmological perturbations are adiabatic. 



turbations occur in single component case 
Then C can be given as 

C = $ + #^. (25) 
P 

It remains constant on super-horizon scales in the adiabatic case even if w is changing. 
Considering the fact that ^ is always equal to — 2$ on super-horizon scales 6l(] and using 
( Fl4|) . we can rewrite Eq. (I25j) as 

3 + 3w — 3a 

Since both w and a axe constants, $ is also a constant. Therefore, we can obtain f e ^ 1 by 
evaluating ($^) and From Eq.([2]), we have the two-point correlation 

_{6E^_ _J_dC/ _ mAT m+1 
{dp) ~ V* ~ V*d(3~ V ' 10 

and the three-point correlation 

3 {8E*) 1 d*U m{m + l)AT^ 
VP ) - ~yr ~ ~ V 2 ' ( 8) 
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Applying the horizon crossing condition k = a/ R = aH and the Poisson equation ([6]), the 
two-point correlation function for is given as 

($1) = {AnG) 2 R\5p 2 k ) = (47rG) 2 Rk~ 3 mAT m+1 , (29) 

while the three-point correlation function is 

($3) = (AttG) 3 R 6 (5pI) = {AnGfk-^m(m + l)AT m+2 . (30) 

Then the non-Gaussianity estimator fx™ 1 can be calculated with the use of Eq. (124jl . It 
turns out 

equi i 5 3 + 3w — 3a (m + 1) 5 3 + 3w — 3a 2w + 1 

Jnl = Y2irGAT m R 2 5 + 3w-3a m = 727rGpR 2 5 + 3w - 3a w + 1 ' ^ ' 

Obviously the amount of non-Gaussianity also depends on both w and a. In the previous 
section, we have concluded that the state parameter should be restricted at w ~ a — 1 or 
w ~ — | in order to obtain a scale invariant spectrum. 

In the case of w ~ a — 1 , we get a small non-Gaussianity 

/X' « 0, (32) 

which is not surprising since it leads to de-Sitter expansion as in the slow-roll inflation 
scenario. 

In the case of w = — |, we have 

/•egmz _ 5 2 — 3a 5 2 — 3a 

/7VL ~ ^vrGpi? 2 4 - 3a ~ 54 4 - 3a ' 1 j 

As a result, if we choose the value of a close to |, a large non-Gaussianity can be obtained. 
Furthermore, we find the non-Gassianity is scale- invariant as well. That is to say, it is not 
suppressed at large scale, which is in contrast to the previous results in string gas cosmology 
and holographic cosmology where the non-Gaussianity is k dependent 
is worthwhile to point out that in this case we have H > 0, which implies a phantom like 
universe. 



24 



25|. Finally, it 



V. CONCLUSION AND DISCUSSION 



In this paper, we have investigated the power spectrum of thermal fluctuations in very 
early stage of viscous cosmology with a special viscous coefficient. Firstly, with the conjecture 

9 



that the origin of the structure of our universe is completely thermal, we have demonstrated 
that a nearly scale invariant spectrum can be achieved when the state parameter and the 
viscous coefficient are properly selected. Different from previous literatures, we have only 
introduced the viscosity of the fluid, without modifying the Einstein equation or the thermo- 
dynamical relations of matter sources. Moreover, the results of thermal scenario under the 
condition of w ~ a — 1 is closely analogous to those obtained in standard inflation scenario. 

Secondly, we analyze the non-Gaussianity of the thermal fluctuations in this scenario. 
It is shown that there will be a suppressed non-Gaussianity in the case w ~ a — 1. It is 
just analogous to the usual inflationary phase where w ~ — 1 is required. However, if we 
choose w ~ — |, a large and k- independent non-Gaussianity can be achieved by tuning a 
near | which implies a phantom like universe. It is of interesting to notice that a phantom 
phase is also required in holographic cosmology and the near Milne cosmology when a large 
non-Gaussianity can be obtained. 

In conclusion, due to the viscosity effects of non-perfect fluids, the no-go result of the ther- 
mal fluctuations can be avoided and a large non-Gaussianity can also be achieved. Therefore 
it is desirable to investigate furthermore the influence of viscosity effects on the structure 
formation and anisotropy of the CMB. 
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